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Abstract 

Possible ways of identification are discussed of an electronic order of higher multipoles such as octupoles and 
hexadecapoles. A particularly powerful method is resonant X-ray scattering (RXS) using quadrupolar resonance 
processes called E2. The characteristic azimuthal angle dependence of Ceo. 7 Lao. 3 Be is interpreted as evidence of 
antiferro-octupole order. For PrRu4Pi2, eightfold pattern against azimuthal angle is predicted if its metal- insulator 
transition is a consequence of a hexadecapole order. In non-resonant superlattice Bragg scattering, hexadecapole 
contribution may also be identified because of absence of quadrupole component. 
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J. Introduction In this paper, we shall discuss how 
to observe electronic orders of higher multipoles such 
as octupoles and hexadecapoles (2 4 -poles). These or- 
ders should of course result in anomaly of specific heat 
as in ordinary phase transition, but otherwise make it 
hard to identify the order parameter. Multipole orders 
have attracted growing interest recently. In particular, 
CexL&i-xTSe exhibits a rich H-T phase diagram with 
phases ranging from I to IV [1,2]. The T^u octupole 
ordering model explains most of the characteristic be- 
haviors in phase IV [3], such as the cusp-like behavior 
in the uniform susceptibility [2], absence of magnetic 
Bragg reflection in neutron scattering [4], huge soft- 
ening of C44 mode in the strain susceptibility [5], and 
minute lattice distortion along [111] [6]. 

On the other hand, PrRu4Pi2 shows a metal- 
insulator transition with breathing type staggered 
lattice distortion [7] . The crystal symmetry remains 
cubic in the low temperature phase [8] . The magnetic 
moment is unlikely to participate in the transition 
which is insensitive to magnetic field, and at which the 
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susceptibility shows no distinct anomaly. The large 
anomaly in the specific heat, on the other hand, re- 
quires involvement of /-electron degrees of freedom. 
Thus it has been pointed out [9,10] that a particular 
type of hexadecapole moment can be a candidate of 
the electronic order parameter in PrRu4Pi2. In order 
to reproduce the CEF splittings below the transi- 
tion, however, additional contribution of hexacontate- 
trapoles (2 6 -poles) should be considered [11]. 

In the following, we shall demonstrate that octupole 
and hexadecapole orders show up as a characteris- 
tic diffraction property of X-ray scattering. Especially 
powerful is identification of a multipole through the 
dependence on azimuthal angle. 

II. Multipoles under CEF We begin with proper- 
ties of multipoles under crystalline electric field (CEF), 
which is specified in terms of the one-body potential. 
Under the tetragonal point-group, which is relevant to 
filled skutterudites, the CEF potential is written as 
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in the standard notation[12] where the superscript c 
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and t mean cubic and tetragonal, respectively. In the 
case of / electrons, their orbital angular momentum I — 
3 gives the upper limit of the rank L — 2 x 3 = 6 for the 
CEF potential. In the cubic point group Oh, the term 
yO\ does not appear. The operator 0¥ is proportional 
to the sum of spherical tensors T^' and T^. The 
particular combinations of with the coefficients Al 
transform as scalars in the point group. These scalars 
appear after the intensity parameter W . 

It is obvious from this form that a scalar in the point 
group is in general a multipole in the spherical sym- 
metry. The anisotropy associated with a multipole is 
probed by such experiments that are sensitive to de- 
tailed shape of wave functions. In microscopic terms, 
the shape specifies distribution of charge and spin of 
electrons. The magnetic octupole moment, for exam- 
ple, consists of distribution of magnetic dipole mo- 
ments that add up to zero around the rare-earth site. 
Both spin and orbital angular momentum can be the 
origin of microscopic magnetic moment. The electric 
quadrupole moment, on the other hand, can be nonzero 
for non-spherical distribution of charges. The hexade- 
capole becomes the lowest multipole moment for cer- 
tain non-spherical distribution of charges. 

In terms of the coordinates r = (x, y, z) on the unit 
sphere for describing angle variables, a set of octupoles 
under the cubic symmetry is given by 

r 4u : x(5x 2 - 3), y{5y 2 - 3), z(5z 2 - 3), (2) 

which transforms in the same way as the set x,y,z. 
Thus the dipole and a part of octupole operators are 
mixed. The irreducible representation is labeled T^. 
On the other hand, another set of octupoles 

r 5u : x(y 2 - z 2 ), y(z 2 - x 2 ), z(y 2 - z 2 ), (3) 

does not mix with dipoles, and are labeled Fs u . Among 
even rank tensors, quadrupole and a part of hexade- 
capole (L = 4) operators mix under the cubic symme- 
try. Namely, a set of hexadecapoles 

r 59 : xy(7z 2 - 1), yz(7x 2 - 1), zx(7y 2 - 1), (4) 

transforms in the same way as xy,yz,zx and belongs 
to the F 5g representation. 

Some cases have been found where the CEF states 
allow these multipoles to fluctuate because of their de- 
generacy, or because of small enough CEF splitting. 
A typical example is the four-fold degenerate CEF 
ground state in CeB6. There is an orbital degeneracy in 
addition to the Kramers degeneracy in this case. Then 
one can expect both magnetic and orbital orderings in 
the ground state, unless quantum fluctuation like the 
Kondo effect blocks the order. Actually a magnetic or- 
der called phase III is realized in the ground state of 
pure CeB6, but the Kondo effect becomes dominant as 
the concentration of Ce decreases in Cei-ajLaajBe. 



III. Non-resonant X-ray scattering from multi- 
poles The most straightforward method to probe a 
multipole order is non-resonant X-ray scattering [13]. 
The scattering amplitude is proportional to the form 
factor 

F(q) = (^exp(iq.r J )), (5) 

3 

where rj is the coordinate of j-th electron in a unit 
cell. The structure factor is nonzero only for Bragg 
vectors or a superlattice vector q = Q involving the 
charge density. For a given q, the plane wave is ex- 
panded in terms of spherical harmonics times spheri- 
cal Bessel functions. Taking the z-axis of spherical har- 
monics along direction of q, we obtain 

F{q) = Y.c n (j 2n { q ))(T^)q, (6) 

n=0 

where c n is a numerical factor, and {Tj^ )q with M — 
and L < 6 indicates that q is taken along the z-axis. 
In the presence of an antiferro-orbital order with Fs g 
symmetry, both T%P and T$ contribute to the scat- 
tering. Furthermore we have introduced the notation 
{jh{q)) = / °° drjL{qr)f(r)r 2 , where f(r) is the radial 
wave function of an / electron, and jl {qr) is a spherical 
Bessel function of order L. The scattering amplitude 
depends on direction of q relative to the cubic crystal 
axis through the spherical tensors. On the other hand, 
the dependence on q = \q\ comes from (jh{q)) Thus 
the scattering intensity gives information about rela- 
tive weight of quadrupole and hexadecapole operators. 
The case of CeB6 has been analyzed in ref. [14] which 
demonstrates that the dominant weight comes in fact 
from the hexadecapole component. 

IV. Resonant X-ray scattering from multipoles 

We now discuss resonant X-ray scattering (RXS). Re- 
cent measurement has obtained information of the or- 
der parameter symmetry in phase IV in Ceo.7Lao.3Bg 
below Tiv = 1.5 K [15]. The superlattice Bragg reflec- 
tions at Q = (3/2, 3/2, 3/2) is taken with the rotation 
axis normal to the [111] surface. Both non-rotated (a- 
a'), and rotated (cr-n') polarization channels have been 
measured at the electric quadrupole (E2) resonance 
near the Ce L2 absorption edge. In the a-a 1 channel, the 
scattering intensity exhibits six-fold oscillation, which 
indicates the occurrence of the electronic order with 
three-fold symmetry along [111]. In contrast, the a-n' 
channel shows three-fold oscillation. 

It is demonstrated in ref. [16] that the Tsu-type oc- 
tupole order can reproduce the azimuthal angle depen- 
dences both in the a-a' and ct-tt' channels. We now 
outline the argument of ref. [16]. The RXS amplitude 
per unit cell is given by the formula: [13] 
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where simplifying approximation has been made that 
all intermediate states (m) have the same energy and 
width F. The energy of the incident photon Tiuj is tuned 
close to the absorption edge A of the relevant atom. 
The amplitude is proportional to 



= (f\e ■ P\m){m\e ■ P\i) 
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where the first term in the right-hand side comes from 
the electric dipole (El) transition, and the second term 
from the quadrupole (E2) transition. Here we have de- 
fined the dipole and the quadrupole operators for each 
ion as P a = e£: t0n Vn and Q a0 = e £* tom r«r£/2, 
and the matrix made from the photon wave and polar- 
ization vectors, X aj j — fc"e /3 /2. The prime represents 
quantities concerning the scattered radiation. The 
maximum rank probed by the E2 transition is L = 4. 
Namely, the triakontadipole (L = 5) and hexacontate- 
trapole cannot be detected by the RXS. 

Table 1 summarizes the azimuthal angle depen- 
dences of multipole scattering in the E2 scattering for 
the case where the crystal coordinate R = (X, Y, Z) is 
the same as the scattering coordinate r. The latter is 
defined by taking the scattering plane as the j/2-plane, 
and the origin of ip along the y axis [16] . Here we have 
used the real tensors for M > defined by 
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The tensor is proportional to the multipole op- 

erator 0|f . If the rotation axis for the azimuthal scan 
is different from the crystal axis Z, one has to relate 
the coordinates r and R by linear transformation. For 
example, Table 2 shows the result for the case where 
rotation axis of the azimuthal scan is [111]. If degen- 
eracy remains in the irreducible representation F, we 
take a linear combination so that the principal axis of 
the multipole is parallel to [111]. 

We define the amplitude /5* 11 ' (a) from the [111] do- 
main for the scattering channel a (= a, 7r). The other 
three equivalent domains fi = [111], [111], [111] have 
different amplitudes fg u (a) in general. The total inten- 
sity is proportional to 



A 
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We take the angle £ = 39° as required from the inci- 
dent photon energy and the scattering wave vector q 



Table 1 

Angle dependence of RXS amplitude of octupolcs T*- 3 -' and 
hcxadecapoles T* 4 ' for the E2 transition. The elevation angle 
is written as£ — tt/2 — with 6 being the incident angle. 
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Table 2 

Angle dependence of RXS intensity with [111] as the rotation 
axis. The order parameter is rank p multipoles with irreducible 
representation T in the Oh group. The origin — is taken 
in the [111]-[112] plane. 
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[15]. We assume that the four equivalent domains are 
equally populated, i.e., w — 1/4. The three- fold oscil- 
lation is explained in terms of domains whose principal 
axis is other than [111]. Each domain gives the scat- 
tering intensity dependent on azimuthal angle. Since 
each domain is rotated by 27r/3 around [111], the sum 
of the intensity aquires the three-fold symmetry inde- 
pendent of details of contribution of each domain. The 
intensity factors are chosen as A aa i — 120, A a7I > = 70. 
The difference between A aa i and A a7I > may be ascribed 
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Fig. 1. Comparison of azimuthal angle dependence of RXS 
between theory and experiment. The symbols are taken from 
Mannix et al. [15]. The dotted and the dashed lines are the 
result for the T5 u -typc octupole order as described in the text. 



using the E2 transition can identify the hexadecapole 
order. On the other hand, the non-resonant scatter- 
ing also has contribution from hexadecapoles as shown 
in eq.(6). Although the monopole contribution with Co 
should dominate the intensity, one may hope to iden- 
tify the contribution from L = 4 since the quadrupole 
component L — 2 should be absent. Identification of 
the L = 4 component provides alternative detection of 
the hexadecapole order. 

In summary, we have discussed that higher multi- 
poles can be detected by RXS in the quadrupole (E2) 
transition. The observed azimuthal angle dependence 
in Ceo.7Lao.3B6 is consistent with the pure octupole 
order with the Es„ symmetry, provided four equivalent 
domains are present in the sample. If our interpretation 
is correct, the three-fold symmetry should not appear 
for a sample with a single domain. A new experiment 
is proposed to probe possible antiferro-hexadecapole 
order in PrRu4Pi2. 
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to extrinsic background from non-resonant contribu- 
tion and/or different energies of intermediate states as 
noted below eq.(7). The fitting with these intensity fac- 
tors reproduces semi-quantitatively the observed data. 
Among all multipoles accessible by the RXS, only the 
E.5u-type octupole yields the maximum at ip — in the 
<t-<t' channel, which is consistent with the observed os- 
cillation [15]. The [111] domain predominates over the 
other domains in the a-a' channel. In the a-n' channel, 
on the other hand, the contribution from the [111] do- 
main is much less than that from the other three equiv- 
alent domains. As a result, a threefold component ap- 
pears with the maximum at ip — 0. Therefore, the Es u - 
type octupole in three-fold axis with the four equiva- 
lent domains explains the RXS result in Ceo.7Lao.3B6. 

V. X-ray scattering from PrRuiP\2 We consider 
how to check possible hexadecapole order transition in 
PrRu4Pi2. According to Table 1, hexadecapoles can be 
probed by E2 scattering. The combination O4 + 50| 
constitutes a scalar in the cubic symmetry. Here O4 
corresponds to T* 4 in Table 1, and give rise to a char- 
acteristic dependence on azimuthal angle. The depen- 
dence cos 4-0 in the a-a' and sin4i/j in the a-n' chan- 
nels are not shared by amplitude for other multipoles. 
In particular, there is no contribution from O4 in the 
a-n' channel. The resulting intensity has the eight-fold 
pattern (oc 1 — cos 8^) with a minimum at ip = 0. The 
a-a' channel will show a four-fold pattern by combi- 
nation with 0®. In this way resonant X-ray scattering 
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